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Coarse-grained model of adsorption of blood plasma proteins onto 
nanoparticles 

Mender Lope^ and Vladimir Lobaskirf^ 

School of Physics, Complex and Adaptive Systems Lab, University College Dublin, Belfield, Dublin 4, 

Ireland 

We present a coarse-grained model for evaluation of interactions of globular proteins with nanoparticles. The 
protein molecules are represented by one bead per aminoacid and the nanoparticle by a homogeneous sphere 
that interacts with the aminoacids via a central force that depends on the nanoparticle size. The proposed 
methodology is used to predict the adsorption energies for six common human blood plasma proteins on 
hydrophobic charged or neutral nanoparticles of different sizes as well as the preferred orientation of the 
molecules upon adsorption. Our approach allows one to rank the proteins by their binding affinity to the 
nanoparticle, which can be used for predicting the composition of the NP-protein corona. The predicted 
ranking is in good agreement with known experimental data for protein adsorption on surfaces. 


I. INTRODUCTION 

When uncoated nanoparticles (NP) enter a living 
organism, they are first exposed to biological fluids, 
which results in a formation of stable or transient NP- 
biomolecule complexes. For large NPs, the biomolecu- 
lar coating is referred to as (protein) corona. It has 
been shown that composition and structure of the corona 
determines the biological reactivity and toxicity of the 
NpPSl 

as well as the NP systemic transport including 
NP uptake into cells. The content of the corona can be 
directly linked to the toxic effects and used to predict 
the toxicity of engineered nanomaterialsIn addition to 
the study of possible hazards, the interest to nanobio in¬ 
teractions is driven by promising applications of NPs in 
food, cosmetics, and medicine.1^1® A quantitative model 
of NP corona formation can facilitate designing of new 
nanomaterials with specific functions. 

The composition of the protein corona and protein 
adsorption kinetics have been studied extensively by a 
broad range of experimental techniques such as fluores¬ 
cence correlation spectroscopy (FCS),^ differential cen¬ 
trifugal sedimentation (DCS) combined with imaging 
techniques,!^ quantita tive l iquid chromatography mass 
spectrometry and dynamic light scatter¬ 

ing (DLS) combined with isothermal titration calorime¬ 
try (ITC)P3I (for a recent review see Ref.!!^. Despite 
the valuable information that experimental techniques 
have provided, there is still much controversy and gaps 
in the physical picture of protein adsorption on NPs: 
disagreement on whether the adsorption is reversible, 
whether proteins change conformation and preserve their 
functionality when complexed with certain particle type, 
whether the corona survives the NP uptake into the cell, 
etc. Undoubtedly, computer simulations can assist ex¬ 
perimental data and reveal molecular scale information 
required to understand the corona formation process. 
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Full atomistic simulation of protein on surfaces have 
already proved useful to advance the understanding of 
molecular interactions that determine the binding of pro¬ 
teins to inorganic nanopart icles.!!^®!] xhe atomistic sim¬ 
ulations are however limited to systems composed of one 
or few proteins and give information well below the time 
scales relevant for the formation of the protein corona. 
A solution to overcome this restriction is to use coarse¬ 
grained (CG) models that reduce the number of interac¬ 
tion sites used in the simulation but keep the required 
molecular information about the proteins and the NPs. 
Some CG models to study the kinetics of the protein 
corona formation have already been proposed (see Vi- 
laseca et al.^^ Bellion et al.^^ Oberle et al^^, as well 
as the section on computer simulations of the review by 
Rabe et alW^, but most of these works use rather sim¬ 
plistic presentation of the proteins and lack molecular 
detail that could be essential for the adsorption kinetics. 

In this work, we develop a CG model that allows us 
to calculate the adsorption energies of arbitrary globular 
proteins onto hydrophobic NPs of arbitrary size. The 
model is built starting from the molecular structure of the 
proteins, and the size of the NP is explicitly included in 
the model. In Section |n| we give a detailed description of 
our model and describe the parametrization process. In 
Section III we show the numerical results on adsorption 
of six most abundant human blood plasma proteins on 
NPs of different radii and charge. Finally, in Section [V| 
we present the conclusions. 


II. MODEL 

A. Coarse-grained protein model 

Our main aim is to design a CG model that would 
reflect NP-protein and protein-protein interactions and 
could be scaled up to simulate multiple biomolecules in 
contact with a NP on relatively long times. Much work 
has been done recently on s yst ematic coarse-graining of 
the proteins (for reviews see^^^^. Based on the previ¬ 
ous experience, we propose a single-bead-per-aminoacid 
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model and consider each protein molecule as a rigid body. 
This model preserves the shapes and sizes of the proteins 
(and therefore, their mobilities and excluded volume ef¬ 
fects) as well as the surface charge distribution, so we 
hope to be able to address their competitive adsorption 
on the NP surface. We use crystal structures of the pro¬ 
teins as obtained from the Protein Data Bank (PDB) 
and place one bead for each aminoacid at the position of 
the corresponding a-carbon atom. Fig. shows an ex¬ 
ample of coarse-graining of protein ai-antitrypsin (AlA) 
as taken from the PDB file ID: 3NE4 and our one-bead- 
per-aminoacid CG model. For the NP, we will consider 
here only spherical homogeneous objects, so that a single 
bead presentation is sufficient. 



FIG. 1. Two representations of the oi-antitrypsin (AlA) 
molecule, (a) A full atoms representation taken from the PDB 
file ID: 3NE4 and (b) the one-bead-per-aminoacid model pro¬ 
posed in this work. 


B. Nanoparticle-aminoacid interactions 

We assume pairwise additivity of all interactions and 
present the net NP-protein interaction energy {U) as a 


sum of the individual interactions of the amnoacids that 
compose the protein with the NP material. Eurthermore, 
U is not only a function of distance from the surface 
of the NP to the center of mass (COM) of the protein, 
dcoM^ but also depends on the protein orientation, which 
is characterized by two Euler angles 0 and 0 (defined in 
the next section). In our model, U{dcoM^ 4^^^) includes 
two contributions: 


N 

U{dcOM, <P,0) = J2 ’ (1) 

i=l 


where N is the number of aminoacids in the protein, 
^Vdw -g Waals interaction of aminoacid i 

with the surface and Uf is the electrostatics interaction 
of the aminoacid i with the surface. 

As we are interested in studying the effect of the size 
of the NP on the adsorption energy, we present the van 
der Waals interaction potential in a form that explic¬ 
itly includes the radius of the NP as a parameter, fol¬ 
lowing the well known Hamaker procedure.^ We start 
from the residue-residue interaction potential proposed 
by Bereau and Deserno,!^ which is based on a modi¬ 
fied 12-6 Lennard-Jones potential. In this model, each 
aminoacid residue is characterized by the hydrophobicity 
index (e^) and we additionally assume that any surface 
segment of the nanomaterial can be modelled in the same 
way as the aminoacid (e^). With these assumptions, the 
interaction between the aminoacid i and a bead of the 
nanomaterial s being at a distance r from each other is 
given by: 


UsAr) = 


4ee 

4ee 

0 , 


+een(l-es,i), r<rc,i, 




rc,i <r < rcut, 
r > rcut, 


( 2 ) 


where Cen is a parameter that scales the interaction en¬ 
ergy, Cs^i is the combined hydrophobicity index of residue 
i and the nanomaterial according to the usual Lorentz- 
Berthelot mixing rules and is given by ag^i 

is the average van der Waals radius of residue i and the 
NP bead, dg^i = {dg -h cri)/2, and rc,i is the position of 
the minimum of the pair potential. We follow the same 
methodology as Ref.^ to define the hydrophobicity in¬ 
dex, which is based on the widely used residue-residue in¬ 


teraction energies proposed by Miyazawa and Jernigan,!^ 
but instead of having a 20 x 20 interaction matrix this 
is reduced to a table of hydrophobicities, one for each 
aminoacid (see Table II irP^. A hydrophobicity index 
0 is assigned to the most hydrophilic residue (LYS) and 
an index 1 to the most hydrophobic one (LEU). At this 
point, we should stress that any other hydrophobicity 
scale can also be used, with the only condition that it 
has to be transformed in such a way that the indexes lay 









3 


between 0 and 1. 

In the above expression, we consider only a small vol¬ 
ume element of nanomaterial, similar to an aminoacid 
in scale, but to further coarse-grain the interaction we 
integrate the energy over a semi-infinite volume of the 


material (with a flat surface) and then over a spherical 
NP. For a flat surface, the interaction potential can be 
expressed in terms of d, the distance between the residue 
center of mass and the closest element of the surface. An 
integration of the 12-6 potential defined in Eq. Q over 
a semi-space gives: 


ur {d) = < 


^esP^s 




2 ) 


(1 






d d^ i^ 

dc,i ^ d ^ dcut? 
d ^ dcut, 


( 3 ) 


where Cqs = ^en, P is the number density of beads in 
the nanomaterial, d is the distance from the residue i to 
the surface, dc,i = (2/5)Although the density p 
seems to be an important parameter scaling the interac¬ 
tion, it is not an independent quantity and therefore is 


not crucial for our method. From fitting the adsorption 
energy to experimental or MD simulation data, we can 
find the composite quantity e^sp (energy density), which 
is sufficient for further calculations. For a nanoparticle of 
radius i?, a similar integration over the particle volume 
gives: 


Up^ir) = { 


4eespcrfi 




(l5r®it:3+63r^it:^+45r=^i?'^+5i?^)cr^ 

(r 2 -i? 2)9 ('^ 2 _^ 2)3 
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(l5r®i?^+63r^i?^+45r2i?'^+5i?®)crf . 15i?^cr® • 


(r 2 -i? 2)9 


(r^-R2)3 


rc,i <r < rcut, 

^ '^cut, 


( 4 ) 


where r is the distance from aminoacid i to the cen¬ 
ter of the nanoparticle. The distance Vc^i corresponds 
to the minimum of the potential and is the 

value of the function {vc, i) as defined in the range 

^c,i ^ < '^cut- We do not show the general expres¬ 

sion for the position of the minimum as it is too bulky. 
The minimum is located at rc,i — R ^ (2/5)at 
R ^ (^s,i and is moved to shorter distances at smaller 
R. The variation, however, is not very large, at i? ^ oc, 
Vc^i — R ^ 0.858374crs^i, ai R = 200(Ts,i it is 0.858375crs^i, 
at R = 20as^i it is 0.858469crs^^, and at R = 2as^i it is 
0.865242crs^i. 

Note that the potential proposed (Eqs. (§ and @) 
will only give a repulsive interaction between a highly 
hydrophilic surface and any aminoacid residue (ie. defin¬ 
ing eg = 0, gives = 0 for all residues). On the other 
hand, assigning a non-zero value for will only change 
the magnitude of the interaction between the aminoacid 
and the surface but not the shape of the potential. In this 
way, the proposed potential is limited to only hydropho¬ 
bic surfaces and cannot reproduce the attraction between 
the beads across the hydration layer. Because of this lim¬ 
itation, we set the value of = 1 for all simulations. As 
an illustration. Fig. shows the proposed van der Waals 
potential for LEU (e^ = 1) and LYS (e^ = 0). Alterna¬ 
tively, a potential that includes hydration effects, such as 


a 12-10-6 Lennard-Jones potential for residue-residue in- 
teraction j^^ l ^^ l or the modified version proposed by Wei 
and KnottP^ to model residue-surface interactions can 
be used to generate a more general interaction poten¬ 
tial. The main drawback of the use of these more refined 
formulas for the potential is that the parametrization is 
more challenging and the applicability of a set of param¬ 
eters can be very narrow. 

The electrostatic interactions in Eq. 0 are modeled by 
placing point charges on the NP surface. This charges 
interact with the charged groups of the protein via a 
Debye-Hiickel potential. The electrostatic interaction en¬ 
ergy between an aminoacid i and all the charges on the 
surface is given by: 

= (5) 

i=l 

where rij is the distance between the residue i and the 
point charge j on the surface, Xb = ^I (47r5o^r^BT) is 
the Bjerrum length, is the Boltzmann constant, T 
the temperature, Cq the dielectric permittivity of vac¬ 
uum, Sr the relative dielectric permittivity of water, 
the charge of residue i, qj the charge of the point charge 
j on the surface, the total number of point charges on 
the surface and Xd the Debye length (defined through 
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FIG. 2. Van der Waals interaction potential divided by CesP 
as a function of the distance of an aminoacid to the adsorbing 
surface with Cs = 1 calculated using Eq. Residue LEU 
{ei = 1) for three NP radii (R = 5 nm, 20 nm and 100 nm) 
and a slab are shown. Residue LYS (e^ = 0) for a NP with 
R = 5 nm. 


(along the negative z-axis). The third DOF is the dis¬ 
tance from the COM to the closest point of the surface, 
dcoM- Here, we sample (j) from 0 to 350° in steps of 10° 
and 0 from 0 to 170° in steps of 10° (note that 0 = 0° 
is equivalent to 0 = 360°, and that 6> = 0° is equivalent 
to 0 = 180°). Instead of obtaining the actual adsorption 
free energy by calculating the potential of mean force for 
all orientations, we only calculate the potential energy U 
(given by Eq. 0), which is the sum of all the pairwise 
interactions between the surface and the aminoacids. As 
the net adsorption energies are expected to be well over 
ksT and, as the proteins are assumed to be rigid, ne¬ 
glecting thermal fluctuations is justified. Note that ref¬ 
erence orientations must be chosen to define the angles 0 
and 0 0°. For the simulations reported in this work, the 
reference orientation of each protein was the PDB con- 
flguratio n use d to build the CG model (more details are 
given in IID). For each configuration (0^, the total 


potential energy is calculated as a function of distance 
of the COM to the surface, U{dcoMi^i^0j). Following 
a similar approach as in Kokh et al^and denoting the 
reaction coordinate dcoM = the mean interaction en¬ 
ergy for any particular orientation is given by: 


= StttjXbCo, with cq is the background electrolyte 
concentration). In practice, the points charges are evenly 
distributed on the spherical surface of the NP using a 
Golden Section spiral algorithm and all points will have 
the same charge qj given by qj = 4:7TaR‘^ where a is 
the surface charge density of the NP and R is the radius 
of the NP. 


C. Orientational sampling and the calculation of the 
adsorption energy 

Here, we are not considering conformational changes 
during the adsorption process and treat the proteins as 
rigid bodies. Although the adsorption process might lead 
to conformational changes, this usually happens at longer 
times than the molecule reorientation on the surface.^ 
Then, the adsorption energies calculated here should pro¬ 
vide a reasonable insight into the kinetics of the NP- 
protein corona formation although may slightly underes¬ 
timate the energy. 

To identify the most favorable orientation of adsorbed 
protein globule (corresponding to the minimum adsorp¬ 
tion energy) we will follow the method suggested by Sun 
et Briefly, a configurational space scan is performed, 
where a systematic rotation of the protein is used to build 
a complete adsorption map. There are three degrees of 
freedom (DOF) that have to be scanned. Fig. shows 
that any point within the protein molecule can be de¬ 
scribed by a position vector from the COM of the pro¬ 
tein. This vector is characterized by two angles: 0 and 0 
and by rotating the molecule an angle —0 about the 2 : di¬ 
rection and then by an angle — ^ -b 180 ° about the y axis 
will make the position vector point towards the surface 


E{(l)i,0j) = -kBT 


X In 




i,ej) Jo ksT J 


(6) 


where a(0^, Oj) is the maximum interaction distance from 
the COM of the protein to the surface for the given ori¬ 
entation. Then the total mean adsorption energy of the 
system, Eadi can be estimated by averaging over all ad¬ 
sorbed states with Boltzmann weightings 

EJ:PijE{cPi,e^) 

i j 

where Pij = exjp[—E{(l)i,0j)/kBT] is the Boltzmann 
weighting factor. 


D. Details of the simulations, parametrization and 
validation 


Due to complexity of blood plasma, here we will only 
consider the molecules that are most likely to affect the 
NP interactions and aggregation and mediate the NP in¬ 
teraction with the cell membranes. The plasma can then 
be modelled a solution of biomolecules in an implicit sol¬ 
vent with a dielectric constant of water and the Debye 
length corresponding to physiological ionic strength, van 
der Waals interactions set to corresponding triplets NP- 
protein-water, or protein-water-protein, and appropriate 
surface charges on the molecules. We selected six rep¬ 
resentative plasma proteins. As mentioned in Sec. |H A[ 
in our CG model each amino acid of a protein is repre¬ 
sented by a single bead located at the a-carbon position. 
The native structures are obtained from the Protein Data 














5 


Residue 

LYS GYU 

ASP 

ASN 

SER ARG 

GLU 

PRO THR GLY 

€i {£) 

0.00 

0.05 

0.06 

0.10 

0.11 

0.13 

0.13 

0.14 

0.16 

0.17 

Oi (nm) 

0.64 

0.59 

0.56 

0.57 

0.52 

0.66 

0.60 

0.56 

0.56 

0.45 

Residue 

HIS 

ALA 

TYR CYS 

TRP 

VAL 

MET 

ILE 

PHE 

LEU 

u {£) 

0.25 

0.26 

0.49 

0.54 

0.64 

0.65 

0.67 

0.84 

0.97 

1.00 

Oi (nm) 

0.61 

0.50 

0.65 

0.55 

0.68 

0.59 

0.62 

0.62 

0.64 

0.62 


TABLE L Normalized hydrophobicities Ci (taken from Table II irP^and ai for each amino acid. The most hydrophilic residue 
has a of 0, while the most hydrophobic has a value of 1. 


Protein PDB ID Abbreviation Weight fraction Molar mass 

in plasma, % in, kDa 


Human Serum Albumin 

1N5U 

HSA 

5.0 

67 

0 ( 1 -antitrypsin 

3NE4 

AlA 

0.24 

51 

02 -macroglobulin 

4ACQ 

A2M 

0.72 

725 

Fibrinogen 

3GHG 

Fib 

0.4 

340 

Transferrin 

2HAV 

Tra 

0.4 

80 

Immunoglobulin G 

3HR5 

IgG 

1.24 

150 


TABLE IL Description of proteins used for the model plasma: the abbreviations used in the text, globule size, and the protein 
abundance in human blood plasma. 



EIG. 3. Definition of the protein orientation, (a) Any “atom” 
of the protein can be described by a position vector from the 
COM, whose orientation is characterized by two angles 0 and 
6. (b) The angles correspond to azimuthal and polar rotations 
(see Sec. |II C] ) that would turn the original vector towards the 
surface (along the negative z-axis). The remaining degree of 
freedom is the distance of the COM to the surface, dcoM- 


Bank, and in Table |TI| we list the proteins under study, 
their PDB IDa from which the CG model were built and 
the abbreviations that will be used in the rest of the text. 
Table El also summarizes their relative content in blood 
and their molar mass. Although these six proteins repre¬ 
sent the most comm on components of the blood plasma, 
recent studie^^^®^ demonstrated that the protein corona 
can include hundreds of different plasma proteins. Fig¬ 
ure shows the CG representations of the six proteins 
chosen for our study. Note that the range of sizes of the 
proteins is very broad, from a big molecule as Fib (about 
10 X 45 nm) to a relatively small molecule AlA (about 


8x4 nm). After the CG model were built from the PDB 
files, the obtained structures were shifted so the COM 
of the molecules was in the origin of the frame of refer¬ 
ence and this structure was defined as the (0 = 0°,^ = 
0°) orientation. With this definition the first residue 
in the sequence of each protein will have the follow¬ 
ing (0,6>) angles: (21.4°,85.2°) for HSA, (101.0°,126.7°) 
for AlA, (193.9°,48.9°) for A2M, (132.1°,46.4°) for Fib, 
(279.6°,140.2°) for Tra and (6.3°,110.8°) for IgG. 



FIG. 4. CG models of the proteins studied in this work. From 
left to right: ai-antitrypsin (AlA), Human Serum Albumin 
(HSA), Transferrin (Tra), Immunoglobulin G (IgG), Fibrino¬ 
gen (Fib) and 02 -niacroglobulin (A2M). The PDB ID struc¬ 
ture from which each CG model was built from are reported 
in Table ini 

All simulation were performed using ESPResSo pack- 
ag^^ and the cutoff for the interaction potential in 
Eq. was set to rent = 0 nm. For all calculations the 
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simulation box was taken big enough to fit the NP and 
the protein. The units of the simulations are: lengths 
(£) in nm, energy (f) in ksT ^ 4.15 x 10“^^ J taking 
a temperature of T = 300 K, unless specified otherwise. 
For the mass unit (A4) we selected the average mass of 
the 20 aminoacids (ca. 110 Da) hence in our simulations 
all aminoacids have a mass of 1. The values of and 
cFi can be found in Table |I] and as mentioned in Sec. HH 
we will only consider hydrophobic NPs with = 1 and 
(Js = 0.35 nm. 


NPs with negative surface charges with charge den¬ 
sity of —0.05 C/m^ as well as neutral NPs were consid¬ 
ered. As explained in Sec. II B[ the charged surfaces are 
modelled by individual point charges. The surface den¬ 
sity of these charged beads {cTc = Ne/Ai^R^) was set to 
nm“^ for all the simulations, which gives e. g. a 
Ne = 100 for a NP of = 5 nm. Then, we assumed 
that each bead carries a charge of — 0.98e, where e is 
the elementary charge. As we are considering physiolog¬ 
ical conditions, we use Xb = 0.73 nm and Xjj = ^ nm. 
Residue charges at these conditions are +e for LYS and 
ARG, —e for ASP and GLU, and -|-0.5e for HIS. The rest 
of the residues are neutral. 

The only free parameter of the model is pCes in Eq. 
and the parametrization was done by systematically 
changing its value to match experimental data of ad¬ 
sorption of Lysozyme on hydrophobic surfaces (octyl- or 
butyl-sepharose) reported by Ghen et alW^. The native 
structure for our GG model of Lysozyme was obtained 
from the PDB ID: 2LYZ. With pees = 1.972/cBT/nm^ we 
obtain a value of —l.^ksT for the adsorption energy (the 
same as the experimental reported value). 

To validate the parametrization, the adsorption energy 
of Myoglobin (PDB ID: IMBN used for the GG model) 
was calculated using the same value of pees obtained from 
the parametrization. In this way, a value of —5.9/c^T was 
found for the adsorption energy of Myoglobin. This value 
is slightly lower that the experimental value of —TSksT 
also reported by Ghen et. but reproduces the trend 
that Myoglobin adsorbs slightly weaker than Lysozyme 
to a hydrophobic surface. 


III. RESULTS 

A. Mean adsorption energies 

Mean adsorption energies for the six proteins calcu¬ 
lated using Eq. 0 as a function of NP radius are shown 
in Eig. Eirstly, the results show that in all cases with 
the exception of A2M the Ead decreases as the radius of 
the NP increases. This trend is due to a combination 
of two factors: (i) increasing R increases the magnitude 
of the van der Waals attraction as shown in Eig. and 
(ii) increasing the radius of the NP also increases the 
surface exposed to the proteins. These effects are more 
pronounced for NPs of R < 100 nm, and after this radius 
the adsorption energies tend to decrease at a smaller rate. 


Eor A2M the curve is non-monotonic and has a minimum 
of the adsorption energy for R = 20 nm. This protein is 
rather big and has a complex structure, which makes the 
effects mentioned above combine in a non-trivial way. 
We see that for all molecules (except Eib) the total Ead 
for the neutral NPs tends to the neutral slab value at 
R = 500 nm showing as expected that for large NPs the 
size has only a small effect on the van der Waals interac¬ 
tions. Eor Eib, more points would be needed for R > 500 
nm to observe how the Ead converges to slab value but 
our results suggest as in the case of A2M that there is a 
minimum in the Eib adsorption energy. 

Secondly, the effect of the electrostatics is smaller in 
magnitude than the van der Waals contribution for the 
surface charge studied here. This can be confirmed by 
noticing that the difference between the Ead of the neg¬ 
atively charged and neutral NPs are significantly smaller 
than Ead itself. In fact, the electrostatic interactions 
modify the adsorption energy by less than SksT per pro¬ 
tein. 

Thirdly, Tra and IgG attract stronger to the nega¬ 
tively charged surfaces, HSA and AlA to neutral sur¬ 
faces while Eib and A2M do not show a clear pattern. 
HSA and AlA are sightly negatively charged (both with 
a total charge of —6e) so overall electrostatic attraction 
dominates over electrostatic repulsion. Eor Tra, the total 
charge is -hl5e but excluding the contribution from HIS 
the total charge is —4e. As the HIS residues contribute 
half a charge, effectively the positive charge of Tra is 
less localized and the molecule behaves as a slightly neg¬ 
ative object. To confirm this observation, calculations 
of the adsorption energy for the same surfaces but with 
the HIS with no charge and with charge +le were done. 
As expected, Tra with uncharged HIS residues attach 
stronger to the a neutral surface, while Tra with pos¬ 
itively charged HIS residues adsorbs stronger onto the 
negatively charged surface. In the case of IgG, the to¬ 
tal charge of the protein is +28e and stays positive even 
without the contribution of the HIS residues (they con¬ 
tribute +12e), so the overall electrostatic attraction is 
greater for the negatively charged surfaces. Eib is slightly 
positive (+3.5e) so it is expected that it attaches stronger 
to the negative surface. This is the case for all radii apart 
from R = 100 nm. Despite that for this radius the Eib 
molecule attaches stronger to the neutral NP, the differ¬ 
ence in the Ead for both surfaces (neutral and charged) 
is only of about 3.5%, which again shows that van der 
Waals contributions dominate over electrostatic interac¬ 
tions. 

Eor A2M which has a total charge of —5.5e, there is 
no clear indication of the charge effect on Ead- As this 
molecule is big, the relative contribution of the charge to 
the adsorption energy is expected to be small compared 
to the van der Waals contribution. Our results agree with 
this prediction as the differences in the Ead are always 
small (less than 5%) with the exception of R = 5 nm. Eor 
this radius, the NP is so small compared to the protein 
that the contact area includes only few aminoacids. In 
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Radius [nm] 

Ranking 

5 

A2M IgG Fib Tra AlA HSA 

20 

Fib A2M IgG Tra AlA HSA 

50 

Fib A2M IgG Tra AlA HSA 

100 

Fib IgG A2M Tra AlA HSA 

500 

Fib IgG A2M Tra AlA HSA 


TABLE III. Ranking of the adsorption energies for the nega¬ 
tively charged surface. For each NP radius, the proteins are 
sorted from left (stronger adsorption) to right (weaker adsorp¬ 
tion) by their value of Ead- 


this case, the electrostatic interactions have larger rela¬ 
tive effect than for the larger sizes. 

Lists of the proteins for each NP size, sorted by the 
adsorption energy, are reported in Table III In all cases, 
the big proteins (Fib, A2M and IgG) are within the three 
those most strongly attached to the NP independently 
of the radius, while the small proteins (AlA, Tra and 
HSA) are in the group of the three molecules with weaker 
adsorption. Also, HSA is the weakest attached protein 
and Fib has always the most negative adsorption energy 
for NPs of R > b nm. 


B. Adsorption energy maps and preferred orientations 

The systematic sampling applied for the calculation of 
the adsorption energies can also be used to identify the 
most favorable orientations for the adsorption. As an ex¬ 
ample, Fig.j^shows a color map of the adsorption energy 
as a function of the angles 0 and 0 for HSA adsorbing 
on a neutral 20 nm-NP. The energy landscape is complex 
in structure showing several connected local minima with 
differences that are less that IksT. This observation sug¬ 
gests that orientational changes in the NP-protein com¬ 
plex after adsorption are likely to occur. Additionally, 
we should note that the map has big areas with adsorp¬ 
tion energies of — or lower, which in practice means 
that more than one orientation gives a relatively strong 
adsorption, so that the proteins will bind to NPs at room 
temperature in various orientations. For the rest of the 
proteins and for the conditions studied in this work, sim¬ 
ilar features are observed in the energy maps. 

To study the effect of the radius of the NP and the 
charge on the adsorption maps, in Fig.j^we show the nor¬ 
malized adsorption energy maps for HSA. The adsorption 
energies have been normalized by rescaling the energies 
such that 0 denotes the the maximum adsorption energy 
while -1 denotes the minimum. This transformation has 
been performed to better illustrate the similarities or dif¬ 
ferences between the different cases. In Fig.[^ each panel 
is for a radius of 5, 20 or 500 nm, respectively, and for a 
neutral or a negative charged surface. A first comparison 
of the structure of the maps in the different panels reveals 
that neither the radius nor the charge density have a ma¬ 
jor impact. A closer inspection of the maps for the same 


charge but different radii (compare Figs. [^,[^ and[^ or 
Fig.[^,[^ and[^) shows that there are only small differ¬ 
ences in the structure of the maps. The largest changes 
are seen for = 5 nm compared to = 20 or 500 nm. 
On the other hand, the effect of the charge on the struc¬ 
ture of the maps is even smaller (compare Fig. with 
or Fig.[^ with[^ or Figs.j^ with[^). We performed 
similar analysis for the other proteins and found that the 
energy surfaces for Tra and IgG are again very stable 
regardless of the surface charge or the radius of the NP. 

A different situation is observed for A2M. Fig|8] shows 
the adsorption energy map of A2M for three radii and 
two charges. The energy maps for = 5 nm show clear 
differences with the maps for R = 20 nm and 500 nm. 
For the smallest radius {R = 5 nm) the maps contain few 
local minima compared to the other two radii. These are 
also more localized than the local minima observed for 
R = 20 and 500 nm. As in the case of the total adsorp¬ 
tion energies, these results can be explained by the size 
and shape of the big A2M molecule. For small NPs, the 
aminoacids that come in contact with the surface in each 
relative orientation are determined by specific patch of 
the protein and as the NPs increases in size these patches 
become larger allowing the NP to interact with a larger 
part of the molecule. The results show that for A2M, a 
NP of R = 5 nm is small enough for the adsorption energy 
to be affected by the local structure of the molecule but 
for a = 20 nm this effect is already lost. With respect 
to the charge, we observe no big difference in the energy 
maps. As mentioned before, the charge has a small effect 
on the total adsorption energy so it is expected that it 
would not dramatically change the energy maps. 


For Fib molecule, we observe a different dependence of 
the energy maps on the NP size. Fig. shows the energy 
landscapes obtained for NPs of radii 5, 20 and 500 nm 
for the neutral and charged surfaces. In this case, dif¬ 
ference can be appreciated from comparing the maps for 
the three radii (compare Fig. [^[^, and[^ or Fig. [^, 
W and[^). Fib is not only a big molecule but it is also 
long, so the size effect explained above for AIM is en¬ 
hanced. The energies for Fib calculated for the other NP 
radius used in this work indicate that for a i? > 50 nm 
the maps do not contain noticeable differences (results 
not shown). To better illustrate this result, in Fig. [lQ|we 
show the most favorable orientations for Fib on a neutral 
surface for four different NP radii {R = 5, 20, 50 and 
100 nm). For the two smallest NPs (Fig. 10r and[T^), 
Fib has its adsorption energy minimum in a configura¬ 
tion where the NP is attached to a tip of the molecule. 
The main difference between the R = b nm and 20 nm 
NP-protein complexes is that the second one is interact¬ 
ing with a larger portion of the tip. Meanwhile, as the 
NP increases in size, Fib tends to adsorb in a sidewise 
orientation (Fig. [1Q| : and [T^). These means that the 
most preferred orientation is the one that corresponds to 
the longest axis of the Fib molecule stretched along the 
surface, which maximizes the number of the aminoacids 
in direct contact with the NP. As for the other proteins, 












FIG. 5. Adsorption energies as a function of the NP radius for the six protein studied and two surface charge densities: 
Negative = —0.05C/m^ and Neutral = no charge. (A) HSA, (B) AlA, (C) A2M, (D) Fib, (E) Tra and (F) IgG. The dashed 
lines show the adsorption energy for the case of a neutral flat surface. 
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FIG. 6. Adsorption energy map for HSA interacting with a 
neutral NP of i? = 20 nm. Energy are is in ksT. 


the surface charge of the NP does not induce noticeable 
differences (compare Fig. with[^ or Fig. with ii 
or Fig. 1^ with|^). For Tra molecule, we find similar 
behavior for the dependence of the energy landscape on 
the NP size to that observed for Fib, with the only dif¬ 
ference that the landscapes do not change anymore for 
R > 20 nm. The effect is again due to the size of the 
molecule, which is not as big as Fib, although elongated 
and thus is distinct from the more spherical AIM or IgG. 


C. Effect of temperature on the adsorption energy 

The results reported until now were obtained at T = 
300 K, which is a temperature commonly used in in vitro 
experiments. Now we also considere the adsorption of 
HSA at a temperature of T = 310 K, which is more 
relevant for in vivo conditions. For the calculations at 
T = 300 K the value of the free parameter of the model 
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FIG. 7. Normalized adsorption energy maps for HSA adsorbing on: (a) Neutral NP oi R — ^ nm. (b) Negatively charged NP 
of i? = 5 nm. (c) Neutral NP of = 20 nm. (d) Negatively charged NP of i? = 20 nm. (e) Neutral NP of = 500 nm. (f) 
Negatively charged NP of = 500 nm. 


{pees = 1.972/cBT/nm^) was obtained from matching the 
Ead calculated using Eg. Q to the value of —I.^UbT 
reported by Chen et as the adsorption energy of 
Lysozyme on hydrophobic surfaces at T = 3 00 K (for 
more details on the parametrization, see Sec. II D| ). In 
the same experimental work, a value of —S.2k bT for ad¬ 
sorption energy of Lysozyme on a hydrophobic surfaces 
at T = 310 K is given which we used to scale the free 
parameter of the model. In this way, for the simulations 
at T = 310 K a value of pCes = 2SEJkBTj nm^ was used. 
Additionally, the parameters for the electrostatic inter¬ 
actions were also changed to account for a temperature 
of 310 K: = 0.72 nm and \b = 0.96 nm. 


Figs. [m and Figs. [m show the adsorption energy 
for HSA as a function of R at the two studied temper¬ 
atures for the neutral and negatively charged surfaces, 
respectively. As expected, the HSA molecule attaches 
stronger at T = 310 K to both surfaces (neutral and 
negative). More importantly, the temperature does not 


change substantially the shapes of the curves, which sug¬ 
gests that the effect of the temperature increase simply 
shifts the adsorption energy toward the more negative 
values. As mentioned in the previous sections, the van 
der Waals interactions dominate over electrostatics and 
as they depend linearly on the parameter ptes the overall 
effect merely reflects this trend. To get a better insight 
on how this energy change depends on the radius of the 
NP, in Fig. EH the ratio between the adsorption energies 
at the two temperatures is shown. For the smaller NP 
{R = b nm) the ratio is lower than for the other cases, in¬ 
dicating that for this radius the effect of the temperature 
on the adsorption energy is more pronounced. For large 
NPs, the ratio is greater and tends to a value of about 
0.88, which is less than the ratio between the energy scal¬ 
ing parameters peek’s for T = 300 K over T = 310 K of 
0.92. Thus, we see that the effect of the temperature 
on the adsorption energies is not only in the uniform de¬ 
crease of the energy, which is not very surprising as Ead 
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FIG. 8. Normalized adsorption energy maps for A2M adsorbing on: (a) Neutral NP of = 5 nm. (b) Negatively charged NP 
of i? = 5 nm. (c) Neutral NP of = 20 nm. (d) Negatively charged NP of i? = 20 nm. (e) Neutral NP of = 500 nm. (f) 
Negatively charged NP of = 500 nm. 


for each orientation defined through an sophisticated in¬ 
tegral (see Eq. (§), so modulating the value of pe^g will 
affect the overall Ead in n non-linear way. As for T = 300 
K, the charge of the surface at 310 K has a small effect 
on the adsorption energies. Finally, we compared the 
adsorption energy landscapes for the two temperatures 
(results not shown) and found that they do not differ 
from each other. 


IV. DISCUSSION 

We can now validate the CG model by comparing 
the predictions with known simulation and experimen¬ 
tal data. Lacerda et. measured the binding asso¬ 

ciation constant for HSA, Fib and a set of y-globulins 
on citrate-coated gold NPs (which can be considered as 
negative moderately hydrophobic NPs). They find that 
for NPs with = 50 nm (this was the biggest in their 


study), the y-globulins are the proteins that adsorb most 
strongly, followed by Fib and finishing with HSA. Our 
rankings reported in Table |III| are slightly different as 
Fib is predicted to attach stronger that IgG. For NPs of 
15 < i? < 30 nm, the experiments show that HSA is still 
the weakest adsorbing molecule, but Fib shows equal or 
largest binding association constant than the q-globulins, 
which agrees with our results. For R = b nm, the experi¬ 
ments show again that HSA has the smallest binding as¬ 
sociation constant while Fib and the 7 -globulins exhibit 
the same affinity. This again fits well with our finding 
as Fib and IgG have a very similar Ead for i? = 5 nm 
{ — IS.SksT for IgG and —IS.SksT for Fib). This discrep¬ 
ancy for the bigger radius are mainly due to that in are 
calculations we are considering only one 7 -globulin while 
the experimental data was collected for a set of proteins 
of similar structure. It is also interesting to compare our 
results with the simulations reported by Vilaseca et al^ 
who studied competitive adsorption of proteins on sur- 
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FIG. 9. Normalized adsorption energy maps for Fib adsorbing on: (a) Neutral NP oi R — ^ nm. (b) Negatively charged NP 
of i? = 5 nm. (c) Neutral NP of = 20 nm. (d) Negatively charged NP of i? = 20 nm. (e) Neutral NP of = 500 nm. (f) 
Negatively charged NP of = 500 nm. 



FIG. 10. Most favorable orientation for the adsorption of Fib on a neutral NP of radius (a) 5, (b) 20, (c) 50 and (d) 100 nm. 
The Fib molecule is color coded the same way in all cases. 


faces. Using CG MD simulations, they found that for a 
flat surface at long times the most abundant protein ad¬ 
sorbed were Fib, then IgG and at last HSA, which also 
agrees with the ranking based on adsorption energies. 
Generally, this phenomenon of adsorbed protein replace¬ 
ment is known in literature as the Vroman effect 


In general, we find bigger proteins adsorb stronger on the 
surfaces, even for small NPs. 

It is important to remark that the adsorption energies 
and rankings calculated in this work can help to predict 
the long-time composition of the NP-protein corona but 
at short times other factors such as the protein sizes, their 
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NP radius (nm) 


FIG. 11. Adsorption energies for HSA as a function of the 
NP radius for two different temperatures, (a) Adsorption 
energies for neutral NPs at T = 300 K (circles) and T — 310 K 
(squares), (b) Adsorption energies for negatively charged NPs 
at T = 300 K (circles) and T = 310 K (squares). The dashed 
lines are the values of Ead for the slab cases at T = 300 K 
(red) and T = 300 K (blue), (c) The ratio 7 = Ead{T = 
300K)/F^ad(T = 310K) for neutral (diamonds) and negatively 
charged (triangles) NPs. 


concentrations and mobilities have to be considered. As 
calculations for any protein with known crystal structure 
can be done easily in our CG model, its results can be 
used as an input for studying competitive adsorption of 
proteins and t he cor ona formation kinetics, such as the 
models in Refs!^^^ 

In Fig. we show a color-coded representation of 
HSA divided in three domains and the most favorable 
orientation for adsorption on a = 500 nm NP (for 
HSA at this radius the NP the adsorption preference is 
equivalent to that on a flat surface). The HSA molecule 
adsorbs in a side-on configuration, in which all three do¬ 
mains are interacting with the surface. Khan et al.^ 
have reported recently a similar result for the adsorption 



FIG. 12 . (a) Domains of HSA molecule, (b) Most favourable 
orientation of HSA on a flat hydrophobic surface. 


of HSA on a hydrophobic surface but using docking sim¬ 
ulations based on full atomistic model. We conclude that 
our CG scheme indeed preserves enough information to 
capture the essential adsorption mechanisms. Therefore, 
our model can be used for a quick search for the pre¬ 
ferred configurations to accelerate docking experiments 
or to study competitive protein adsorption in plasma or 
other protein solutions.^ 

The main limitation of our model is that as we consider 
the proteins to be rigid bodies, conformational changes 
are not allowed and in some cases thi s can be an impor¬ 
tant factor for the adsorption processThis assump¬ 
tion can be relaxed by e.g. using elastic network derived 
from the principal component analysis of the molecule 
dynamicJ^ or a Go-Type model (se^^ for a review on 
CG models of proteins). Another limitation is the in¬ 
ability to address properly the interaction of hydrated 
residues/surfaces. This deficiency can be corrected by 
introduction of more accurate potentials of mean force 
that take into account the water and surface structu ring, 
which can be derived from all-atom MD simulations 
We are currently working on combining the PMFs from 
atomistic MD simulations with our model of proteins. 


V. CONCLUSIONS 

In this work, we presented a CG model for calcula¬ 
tion of the adsorption energies of globular proteins on 
hydrophobic charged NPs. The proposed method was pa¬ 
rameterized and validated against more detailed simula¬ 
tion and experiments. The model can be applied for eval¬ 
uation of binding energies for arbitrary plasma, cytoso¬ 
lic or membrane proteins with known structure, ranking 
them by binding affinity to the NP and predicting the 
content of NP protein corona. We performed a study 
of adsorption of six common blood plasma proteins onto 
hydrophobic NPs. We found that the NP surface charge 
has a small effect on the adsorption energies in compar¬ 
ison to van der Waals interactions between the residues 
and the surface. We also found that the charge of the 
NP does not have a major influence on the orientation, 
in which the proteins are to be mostly likely adsorbed. 
On the other hand, we showed that the size of the NP 
has a pronounced effect on the adsorption energy maps. 
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as the curvature of the NP determine the sections of the 
protein that can come in contact with the NP surface. Fi¬ 
nally, we note that, as the methodology presented here is 
computationally efficient and gives consistent predictions 
for the structure of NP-protein complexes, it can be used 
as a part of multiscale modelling methodologies or a sup¬ 
plement to more sophisticated approaches to modelling 
bionano interactions. 


ACKNOWLEDGMENTS 

The presented research has been funded by EU 
FP7 collaborative grant, project 310465 (Membrane- 
NanoPart). 

^M. Monopoli, C. Aberg, A. Salvati, and K. A. Dawson, Nat. 
Nanotechnol. 7 , 779 (2012) 

Lynch, K. A. Dawson, and S. Linse, Sci. Signal. 2006 , pel4 
(2006) 

^T. Cedervall, I. Lynch, S. Lindman, T. Berggard, E. Thulin, 
H. Nilsson, K. A. Dawson, and S. Linse, Proc. Natl. Acad. Sci. 

U. S.A. 104 , 2050 (2007). 

"^S. Lindman, 1. Lynch, E. Thulin, H. Nilsson, K. A. Dawson, and 
S. Linse, Nano Lett. 7 , 914 (2007). 

^L. Allen, M. Tosetto, 1. Miller, D. O’Connor, S. Penney, 1. Lynch, 
A. Keenan, S. Pennington, K. Dawson, and W. Gallagher, Bio¬ 
materials 27 , 3096 (2006). 

®P. Kamath, A. Fernandez, F. Giralt, and R. Rallo, Curr. Top. 
Med. Chem. 15 , 1930 (2015). 

^O. V. Salata, J. Nanobiotechnol. 2 , 3 (2004) 

®M. Rahman, M. Ahmad, 1. Kazmi, S. Akhter, M. Afzal, 
G. Gupta, and S. V.R., Curr Drug Discov Technol. 9 , 319 (2012). 
9Z. Wang, G. Niu, and X. Chen, Pharm Res. 31 , 1358 (2014) 
^°C. Roecker, M. Poetzl, F. Zhang, W. J. Parak, and G. U. Nien- 
haus. Nature Nanotech. 4 , 577 (2009) 

^^P. M. Kelly, C. Aberg, E. Polo, A. O’Connell, J. Cookman, J. Fal¬ 
lon, Z. Krpetic, and K. A. Dawson, Nature Nanotech. 10 , 472 
(2015) 

C. Silva, M. V. Gorenstein, G.-Z. Li, J. P. C. Vissers, and 
S. J. Geromanos, Mol. Cell. Proteomics 5, 144 (2006) 

^^S. Ritz, S. Schoettler, N. Kotman, G. Baier, A. Musyanovych, 

J. Kuharev, K. Landfester, H. Schild, O. Jahn, S. Tenzer, and 

V. Mailaender, Biomacromolecules 16 , 1311 (2015) 

^"^S. Winzen, S. Schoettler, G. Baier, C. Rosenauer, V. Mailaender, 

K. Landfester, and K. Mohr, Nanoscale 7 , 2992 (2015) 

^^P. del Pino, B. Pelaz, Q. Zhang, P. Maffre, G. U. Nienhaus, and 

W. J. Parak, Mater. Horiz. 1, 301 (2014) 


^®G. Brancolini, D. B. Kokh, L. Calzolai, R. Wade, and S. Corni, 
ACS Nano 6, 9863 (2012) 

^^F. Ding, S. Radic, R. Chen, P. Chen, N. Geitner, J. Brown, and 
P. Ke, Nanoscale 5, 9162 (2013) 

^®S. Khan, A. Gupta, and C. Nandi, J. Phys. Chem. Lett. 4 , 3747 
(2013) 

^^F. Tavanti, A. Pedone, and M. C. Menziani, New J. Chem. 39, 
2474 (2015) 

^^P. Vilaseca, K. Dawson, and G. Franzese, Soft Matter 9 , 6978 
(2013) 

^^M. Bellion, L. Santen, H. Mantz, H. Hoehl, A. Quinn, A. Nagel, 
C. Gilow, C. Weitenberg, Y. Schmitt, and K. Jacobs, J. Phys. 
Condens. Matter 20, 404226 (2008) 

^^M. Oberle, C. Yigit, S. Angioletti-Uberti, J. Dzubiella, and 
M. Ballauff, J. Phys. Chem. B 119 , 3250 (2015) 

^^M. Rabe, D. Verdes, and S. Seeger, Adv. Colloid Interface Sci. 
162 , 87 (2011) 

24V Tozzini, Curr. Opin. Struct. Biol. 15 , 144 (2005) 

^^S. Takada, Curr. Opin. Struct. Biol. 22 , 130 (2012) 

26w. G. Noid, J. Chem. Phys. 139, 090901 (2013) 

C. Hamaker, Physica 4, 1058 (1937) 

28t. Bereau and M. Deserno, J. Chem. Phys. 130, 235106 (2009) 
^®S. Miyazawa and R. Jernigan, J. Mol. Biol. 256, 623 (1996) 

^°Y. Kim, C. Tang, G. Clore, and G. Hummer, Proc. Natl. Acad. 
Sci. USA 105 , 12855 (2008) 

Kim and G. Hummer, J. Mol. Biol. 375 , 1416 (2008) 

^^S. Wei and T. Knotts, J. Chem. Phys. 139 , 095102 (2013) 

^^M. Agashe, V. Rant, S. Stuart, and R. Latour, Langmuir 21 , 
1103 (2005). 

^^Y. Sun, W. Welsh, and R. Latour, Langmuir 21, 5616 (2005) 
^^D. Kokh, S. Corni, P. Winn, M. Hoefling, K. Gottschalk, and 
R. Wade, J. Chem. Theory Comput. 6, 1753 (2010) 

^®A. Lesniak, A. Campbell, M. P. Monopoli, 1. Lynch, A. Salvati, 
and K. A. Dawson, Biomaterials 31, 9511 (2010). 

^^H. Limbach, A. Arnold, B. Mann, and C. Holm, Comput. Phys. 
Commun. 174 , 704 (2006) 

^^W. Chen, H. Huang, C. Lin, F. Lin, and Y. Chan, Langmuir 19 , 
9395 (2003) 

^^S. Lacerda, J. Park, C. Meuse, D. Pristinski, M. Becker, 
A. Karim, and J. Douglas, ACS Nano 4 , 365 (2010) 

^°L. Vroman, Nature 196, 476 (1962) 

'^^C. A. LeDuc, L. Vroman, and E. F. Leonard, Ind. Eng. Chem. 
Res. 34 , 3488 (1995). 

Ortega-Vinuesa and R. Hidalgo-Alvarez, Biotechnol. Bioeng. 
47 , 633 (1995). 

"^^M. Holmberg and X. Hon, Langmuir 25 , 2081 (2009) 

‘^‘^S. Koehler, F. Schmid, and G. Settanni, NIC Series 47 , 117 
(2014). 

"^^E. Brandt and A. P. Lyubartsev, Biophys. J. 106, 208a (2014). 
"^^E. Brandt and A. P. Lyubartsev, J. Phys. Chem. C (2015), 
10.1021/acs.jpcc.5b02670 


